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Abstract —Millimeter wave (mmW) bands between 30 and 
300 GHz have attracted considerable attention for next- 
generation cellular networks due to vast quantities of available 
spectrum and the possibility of very high-dimensional antenna ar¬ 
rays. However, a key issue in these systems is range: mmW signals 
are extremely vulnerable to shadowing and poor high-frequency 
propagation. Multi-hop relaying is therefore a natural technology 
for such systems to improve cell range and cell edge rates 
without the addition of wired access points. This paper studies the 
problem of scheduling for a simple infrastructure cellular relay 
system where communication between wired base stations and 
User Equipment follow a hierarchical tree structure through fixed 
relay nodes. Such a systems builds naturally on existing cellular 
mmW backhaul by adding mmW in the access links. A key 
feature of the proposed system is that TDD duplexing selections 
can be made on a link-by-link basis due to directional isolation 
from other links. We devise an efficient, greedy algorithm 
for centralized scheduling that maximizes network utility by 
jointly optimizing the duplexing schedule and resources allocation 
for dense, relay-enhanced OFDMA/TDD mmW networks. The 
proposed algorithm can dynamically adapt to loading, channel 
conditions and traffic demands. Significant throughput gains 
and improved resource utilization offered by our algorithm over 
the static, globally-synchronized TDD patterns are demonstrated 
through simulations based on empirically-derived channel models 
at 28 GHz. 

I. Introduction 

Millimeter wave (mmW) networks are an attractive can¬ 
didate for beyond 4G and 5G cellular system evolution. 
Such systems can potentially offer tremendous increases in 
bandwidth along with further gains from highly directional 
antenna arrays CD, 0. However, a key issue in these systems 
is cell coverage and range due to the extreme susceptibility of 
these high-frequency signals to shadowing and high isotropic 
propagation loss. Given these range limitations, multi-hop 
relaying is a natural technology in the mmW space o, m. 
Furthermore, multi-hop relaying is particularly attractive for 
mmW systems since many cellular systems already use mmW 
backhaul and thus the addition of access links will naturally 
give rise to multi-hop systems. 

Most mmW systems designs have assumed a time-division 
duplex (TDD) structure to fully exploit beamforming and 
eliminate the need for paired bands. In current TDD cellular 
standards, such as TD-LTE, all subframes are globally syn¬ 
chronized with base stations transmitting in one common set 
of DL time slots and User Equipment (UEs) in transmitting 
in the complementary UL set 0. The DL/UL transmission 
mode patterns are essentially static and cannot be adjusted for 
load balancing or changing channel conditions experienced 


by mobile nodes. Additionally, in the current LTE Release 
10 specifications El, in-band relay communication between 
the “donor” eNodeB and Relay Node (RN) can occur only 
in designated subframes, which, as we will show in Sec¬ 
tion |V| could result in the wireless backhaul being severely 
bottlenecked. In the case of such traditional relay-enhanced 
networks, specific subframes must be reserved for the purpose 
of Inter-cell Interference Coordination (ICIC): Since eNBs 
typically transmit at a much greater power than the RNs, 
if these nodes transmit together, the RN signals may be 
overwhelmed by interference from the donor cell. 

Such static, synchronized duplexing may not be necessary 
and, in fact, may be particularly disadvantageous for mmW 
systems and wireless systems that use high-gain, directional 
antennas. In these systems, interference from transmitters can 
be isolated even if there are significant power disparities (as 
found in Q). In this work, we thus consider a dynamic 
duplexing scheme where individual links can select their own 
transmit-receive duplexing pattern. Specifically, we consider 
a system where subframes are synchronized network-wide in 
time, but the transmit/receive selections in each subframe can 
be made on a link-by-link basis - a feature uniquely available 
in the mmW range due to directional isolation. This flexibility 
enables the duplexing pattern to be dynamically optimized 
to current traffic loading and channel conditions. In addition, 
the duplexing can be adapted to local topological constraints. 
This adaptation is particularly valuable since the number of 
hops and their capacity are likely to vary significantly due to 
different cell sizes, propagation obstacles and availability and 
quality of wired backhaul. 

Related work: There is now a large body of work in 
optimization, scheduling, power control and relay selection 
in OFDMA/TDD cellular networks ll8l l9l. iflQl. ifTH. fl2l. 
ca, including several that propose dynamic TDD algorithms 
designed to take advantage of the new LTE-B enhanced Inter¬ 
ference Mitigation and Traffic Adaptation (elMTA) capabili¬ 
ties ifT4i . ifBi . m. These works focus mainly on scheduling 
and ICIC in the context of 4G microwave networks under 
the assumption of an interference-limited regime, whereas we 
assume constant interference due to directional isolation. We 
also take the approach of centralized scheduling, which is 
in contrast to works on distributed MAC schemes for mesh 
networks, as in fTT) . Also, while our analysis is based on 
simulation, we point out that stochastic geometry analysis of 
self-backhaul mmW networks has recently appeared in ms 
as well as a scaling law analysis in CD. 


II. System Model 



Figure 1: Directional multi-hop cellular network with millime¬ 
ter wave backhaul and access links 

A. Network Topology 

We consider an OFDMA-TDD cellular network with di¬ 
rectional smart antennas and multi-hop, in-band relaying as 
shown in Fig. [T] We assume a tree-structured network, where 
an eNodeB base station (denoted BS ) with wireline back¬ 
haul provides a root from which connections go to Njje 
user equipment mobiles (UEi, ..., UEn ue ) via Nrn relay 
nodes RNi ... RNjy RN that are self-backhauled, decode-and- 
forward stations. RNs are essentially indistinguishable in oper¬ 
ation from the BS but utilize wireless backhaul to the primary 
BS, which, unlike the BS’s wired backhaul connection, is sub¬ 
ject to the impairments of the mmW channel. For exposition, 
we assume a two-hop network, although all the methods apply 
to multihop networks as well. Each user is associated with 
either the BS (i.e. direct-link) or a relay RNj. The rate on the 
link between RNj and its parent BS is denoted Ri 1:j and Ri jtl 
for the downlink (DL) and uplink (UL), respectively, The rates 
on the DL and UL links between UEi and its parent BS or 
RN are likewise written Ri j i and Ri i}j - 

We define T to be the set of flows in the network. Users 
have one DL flow and one UL flow, which gives Njr = 2Njje 
total flows in E, but it can be easily extended to allow for 
multiple differentiated types of traffic. Flows have a throughput 
Rf and utility Uf(Rf). We assume a standard proportional- 
fair metric l20l 

U f (R f ) = \og(R f ), (1) 

although other concave utilities may also be used. 

Users are scheduled in a series of epochs or frames of 
period Tf, which are further subdivided into N s f subframes 
of period Tsf = Tf/N s f. From the perspective of each node, 
each subframe can be designated for DL or UL transmission 
indicating whether the transmissions in that subframe are away 
from or toward the BS root of the tree. A subframe can also 
be muted (i.e. unutilized). Within each subframe, we assume 


Orthogonal Frequency Division Multiple Access (OFDMA) is 
employed, allowing multiple users to be allocated orthogonal 
frequency resources within same subframe. Specific OFDMA 
parameters, discussed in Section [V| are derived from the LTE- 
like mmW system design proposed in ED. 

In contrast to the semi-static and globally synchronized 
TDD configurations supported by relay-enhanced TD-LTE 
networks, we allow each individual BS and relay node to 
dynamically select the transmission mode of each subframe. 
The assignment is centrally-coordinated through control mes¬ 
saging from the BS, however the operation of the specific 
MAC protocol is beyond the scope of this work. 

Additionally, we make the following assumptions about 
resource assignment and the network, in general: 

i. Half-duplex - Transmission and reception cannot occur 
simultaneously (i.e. during the same subframe) between 
pairs of adjacent nodes. 

ii. Constant-interference - Wireless links behave nearly like 
point-to-point links due to the spatial isolation of direc¬ 
tional beams. This is the key insight that allows us to 
analyze each routing tree, rooted at the BS, separately, 
and schedule resources without need for interference 
coordination. Based on the findings in 0, we assume that 
high-gain, directional beamforming, enabled by multi¬ 
element antenna arrays, provides for minimal interference. 
However, receiving nodes may experience some small 
but non-negligible interference power. To be conservative, 
the interference at each node is taken to be the power 
assuming all nodes in the network were to be simulta¬ 
neously transmitting at full power, with averaged TX/RX 
beamforming directions. 

iii. Single-stream - A node can transmit to only one re¬ 
ceiver in the same time slot and frequency resources. In 
later work, we shall consider BS and RN nodes capable 
of multi-user Space Division Multiple Access (SDMA), 
which can more efficiently utilize the degrees of freedom 
of the channel than standard OFDMA j22l . 

iv. Constant-channel - As in m, we assume that shadowing 
and other large-scale channel parameters vary slower than 
the duration being simulated. Since resource allocation 
is only modeled over a small sequence of frames, these 
stochastic parameters remain constant for the entire se¬ 
quence. Small-scale fading is accounted for in our formu¬ 
lation of achievable rate. 

v. Single-path and static routing - Each RN has a single 
link back to the BS and UEs are associated only with the 
BS or RN that offers the least path loss. The topology of 
the network is therefore fixed over the entire sequence 
of frames. Furthermore, it is assumed that relays are 
deployed by the operator as to guarantee a minimum 
capacity for the mmW backhaul link between the RN and 
wired BS, meaning that UEs will not associate an RN that 
cannot provide a quality backhaul connection. 





B. Channel Model 


For the simulation, we adopt the distance-based path loss 
model for a dense, urban environment at 28 GHz developed 
in m, which itself is empirically derived from extensive 
measurements taken in New York City (231, EH- Each link 
from an access point (BS or RN) to a UE are randomly in 
one of three states: non-line-of-sight (NLOS), line-of-sight 
(LOS) or outage (i.e. the signal is blocked completely by some 
obstruction in the environment). For nodes with a NLOS path 
to their selected transmitter, i.e. the serving BS or RN, the 
path loss is represented by ([2]), 

PL dB {d) = a + /310log 10(d) + £ £ ~ V(0, a 2 ), (2) 

where d is in meters and the parameters a , (3 and a are 
provided in table [T] - see £9 for details. 

As earlier noted, UE nodes initially select their serving 
BS or RN node with the least path loss. The long-term 
beamforming gain is then applied, which is a function of the 
expected gain given the optimal TX and RX beamforming 
vectors, which, in turn, depends on the number of antenna 
elements (jitx and urx) in the A/2 planar antenna array. 
We again defer to (7) for further details on how the long-term 
beamforming gain is calculated, as this subject is outside the 
scope of this paper. 


C. SINR and Rate Calculation 


As discussed in the Introduction, for the purpose of the 
optimization, we assume a conservative interference model 
where the Signal-to-Interference-Noise Ratio (SINR) from any 
transmitting node i to a receiving node j is computed as 


SINR^- 


PrX,i 


int,max,i 


N,i 


(3) 


where PRx,ij is the received power from node i to node j 
assuming the path loss and optimal long-term beamforming 
gain; Pint,max,i is the maximum average interference and Px,i 
is the thermal noise. The maximum average interference is 
computed assuming that (a) all other nodes are transmitting 
and (b) each of the other nodes select a beamforming direction 
to a random node that it could possibly transmit to. When a 
node such as a BS or RN transmits to multiple other nodes, the 
transmit power is scaled by the bandwidth allocated to each 
of the receivers. Given the SINR and bandwidth allocation, 
we assume that maximum data rate for (DL or UL) link lj is 
calculated based on the formulation in [25], which abstracts 
the achievable capacity for LTE/LTE-A systems as 


Cij = vWij log 2 (1 + io 0 - 1 (SINR i j-A 1 „..)) i (4) 


where Wij is the bandwidth allocated to the link, r] is 
the bandwidth overhead and A i oss accounts from loss from 
Shannon capacity due to small-scale fading, channel coding 
and inaccuracy of Channel State Information (CSI). Following 
171 , we take r] = 0.8 and A i oss = 3 dB. 


III. Optimization Formulation 
We consider the problem of maximizing the sum utility, 

U(R):=J2Uf(Rf ), ( 5 ) 

f 

where R is the vector of the flow rates Rf , / E T . We 
divide the frame into N s f subframes, and in each subframe, 
t = 1,..., N s f, and each link £, we let x\ E {0,1} be the 
binary variable indicating whether there is a transmission on 
that link in that subframe. Determining the binary variables x \, 
in effect, determines the TDD scheduling pattern. To enforce 
the half-duplex constraint, we require that in each subframe t , 

x\ + x\, < 1, W = (p, n), £' = (n,c), (6) 

x\ + x\, < 1, Ml — (n,p), £' = (c,n), (7) 

for all links i = (p, n) and (n,p) between a node n and its 

parent p and £' = (n, c ) and (c, n), which the links from n to 

child c. In addition to the transmission schedule variables, the 
optimization will also allocate bandwidths W\ for each link £ 
in subframe t with a maximum bandwidth constraint 

X (8) 

le£(n) 

where Wmax is the total available bandwidth and C(ri) are the 
set of links transmitted from node n. We additionally require 
that W\ = 0 when x\ — 0. Given the bandwidth allocations 
W\, determines the link capacity C\ from ([?]). We let R^^ 
be the rate allocated to flow / on link £ in subframe t which 
must satisfy the flow and capacity constraints 

Kf,i — Ch ( 9 ) 

Rf<EEt=\Rh- do) 

Observe that we have assumed that an arbitrary bandwidth 
fraction can be allocated. Typical OFDMA systems such as 
LTE only permit scheduling of frequency-domain resources 

with the granularity of discrete subcarriers or groups of 

subcarriers (i.e. Resource Blocks). 

IV. Proposed Algorithm 

The problem as formulated is a binary mixed-integer opti¬ 
mization, which is non-convex. A number of algorithms and 
heuristics have been proposed for OFDMA-TDD scheduling, 
many of which employ mixed-integer methods to optimize 
over the integral search space formed by the subcarrier and 
time slot indices rm da, da. These techniques are de¬ 
signed to deal with the computationally hard problem posed 
by an inteference-limited network. As shall be demonstrated 
in the following section, two key aspects of our problem allow 
us to reduce the size of integral search space. Firstly, subframe 
allocation can be performed individually for each BS and RN 
without regard to interference at other nodes. Secondly, the 
hierarchical structure of the network is exploited to identify 
bottlenecks and strategically reassign resources to links to 
iteratively improve utility. 



The suboptimal, greedy algorithm in [T] operates by re¬ 
cursively updating slot configurations on each subtree. The 
procedure TestSubtree takes as arguments, a root node n, 
an initial TX/RX mode indicator matrix X = (x\) for all 
nodes i in the subtree rooted at n and an optimal utility 
for that allocation. Note that x\ represents whether a node 
can be transmitting (x\ = +1) or receiving (x\ = —1) in 
t, which differs from the notation introduced earlier for the 
link allocation xi. In this scheme, the link allocation xi is 
determined by whether one adjacent node is transmitting while 
the other is receiving. 

x\ = l ,W=(p,n) x t p = +1 Ax t n = -1, (11) 

x\ = l ,W=(n,_p) Xp = -1 A x f n = +1, (12) 

The procedure returns a tuple (X, U) of the improved 
duplexing schedule and utility for that subtree. The procedure 
operates recursively, where at each level of the tree, it attempts 
to allocate additional subframes in either the DL and UL (by 
switching the TX/RX mode of nodes at that level) and calling 
the inner optimizer, here denoted as Opt(), to test if utility has 
improved. The inner optimizer utilizes the fact that, for any 
candidate duplexing schedule, the optimization problem is the 
maximization of a concave function with linear constraints. 

We initially call TestSubtree() with the root node n = BS 
to perform the following breadth-first procedure on the entire 
tree. In lines 3-10, we make a recursive call to reschedule each 
RN subtree with the initial X. This tests the case where the 
current configuration is optimal for links adjacent to the BS but 
possibly not for subtrees rooted at a RN. We then reschedule a 
UL subframe in the DL (i.e. x new = +1) at the BS by calling 
Realloc in line 13. We naturally set each UE in the set of 
children C(n) of the root node (the BS, in this case) to RX 
mode in the same t; otherwise, these links would be unutilized 
in this subframe. When Realloc is called on a RN, it selects 
the first t for rescheduling such that its transmission mode is 
the same as its parent. If no such t exists, it simply selects 
the next subframe in succession. The inner optimizer is then 
called, which returns the intermediate utility U ts t- While U ts t 
may be less than the initial utility U at this point, there is still 
the opportunity for it to improve by rescheduling subframes 
in each relay subtree (lines 18-29) (therefore, the procedure 
inside the while loop continues as long as the intermediate 
utility Utst is increasing). In the base case, where the base of 
the tree has been reached or we have returned from testing all 
subtrees, we then attempt to improve utility by rescheduling 
slots in the UL using the same recursive procedure. 

V. Simulation Methodology & Results 

Apart from the parameters taken from 0, we base our 
simulations on the 3 GPP urban microcell (UMi) model and 
guidelines laid down by the 3 GPP and ITU-R for modeling 
LTE/LTE-A networks [26], [27 ]. The use of 1ms subframes is 
based on the Samsung mmW design, presented in l2ll . 

We design our simulation scenarios to compare the per¬ 
formance of our dynamic TDD algorithm with the static TDD 
and relay configurations supported by LTE. For the static case, 


Algorithm 1 Dynamic TDD scheduler 
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procedure TestSubtree( n, X, U) 

Utst CO 

for r G C(ri)\r G TZAf do 

(Ul^Xut) <- TESTSUBTREE(r, X, U) 

if U[ st > U then 
Utst <- U[ st 
u <- Utst 

X <- Xtst 

end if 
end for 

for x new G {+1, -1} do 
while Utst > U do 

X tst ^ REALLOC(X,?7<,X n e^/;) 
for {u G C(n)\u G U£} do 

X tst ^ Realloc(X£ S £, , u, x new ) 
end for 

Utst Opt (X t st) 

for {r G C(n)\r G 1ZM} do 

X tst ^ REALLOC(X,'r, Xnew') 

(Ulrt, Xtst) TESTSUBTREE(r, X tat , U) 
if Ul st > Utst then 
Utst <- U[ st 
Xst <- Xtst 
end if 
end for 

if Utst > U then 
U <r- U t st 
x 4- Xtst 

end if 
end while 
end for 
return (X, U) 
end procedure 


we test every valid combination of the TDD patterns in table 
4.2-2 of f28l and the possible reserved eNB-RN transmission 
subframes in table 5.2-2 of ea. For each static configuration, 
we run the inner OFDMA optimization to provide a fair 
comparison with the dynamic scheduler. 

We test a random network topology with 2 and 4 relay 
nodes. For each scenario and topology, a random drop is 
performed where nodes are placed uniformly randomly in the 
area, with the constraint that RNs must be over 50m from the 
wired BS. Note also that for RN-to-BS links, we assume that 
the deployment of RNs is planned by the operator to ensure 
a LOS backhaul connection. 

The marked gains offered by optimized DTDD scheduling 
over static LTE TDD relaying are shown in Figures [3a] and 
|3b[ Also from table [HI] we see how DTDD provides gains of 
over 1.5x in average downlink rate for the 2-relay and 4-relay 
cases, respectively. We find an even more distinct improvement 
in average uplink flow rate: 3x and 4.5x for the 2- and 4- 
RN scenarios. Improvement in median rates are dramatic, 








Table I: Network model parameters 


i 



Parameter 

Value 

Antenna/RF 

TX power 

P TX = 30 dBm (BS and RN) 

P TX = 20 dBm (UE) 

Antenna 

riant = 64, 8x8 A/2 planar array (BS and 
RN) 

riant = 16, 4x4 A/2 planar array (UE) 

Noise figure 

NF tx = 5 (BS and RN) 

NF tx = 7 (UE) 

Channel 

Bandwidth 

W = 1000 MHz 

Carrier freq. 

f c = 28 GHz 

Path loss 

a = 72.0, P = 2.92, a = 8.7 dB (NLOS) 
a = 61.4, 0 = 2, a = 5.8 dB (LOS) Q 

LOS-NLOS- 
outage prob. 

aout = 0.0334m -1 , b ou t = 5.2, a; os = 
0.0149m -1 0 

TDD 

Frame period 

Tf = 10 ms 

Subframe period 

T s = 1ms 

Subframe 

allocation 

Case 1: Static with reserved BS ++ RN SFs 
Case 2: Dynamic 

General 

Area 

400m 2 (wrap-around distance calculation) 

Inter-site 

Distance (ISD) 

RNs placed at > 50m from BS 

Number of nodes 

Case 1: n^s — 1? rijiN — 2 ,n\jE — 10 
Case 2: ubs — 1? riBN — 4 ,njjE — 10 

Traffic 

Full buffer 

Number of drops 

A drop = 25 


Table II: Mean fractions of LOS, NLOS, and outage UEs 



LOS 

NLOS 

Outage 

2RN 10UE 

0.388 

0.600 

0.012 

4RN 10UE 

0.536 

0.460 

0.004 


Table III: UE flow mean, median and cell edge (worst 5%) 
rate for dynamic and static TDD (units in Mbps) 



Mean 

Median 

Cell edge 

DL 

UL 

DL 

UL 

DL 

UL 

2RN 

DTDD 

31.73 

26.65 

20.86 

10.76 

0.35 

0.19 

static 

21.20 

8.80 

7.70 

0.63 

0.23 

0.15 

gain 

1.50x 

3.03x 

2.71x 

17.24x 

1.57x 

1.28x 

4RN 

DTDD 

48.58 

51.04 

39.45 

28.27 

0.68 

0.37 

static 

29.61 

11.32 

8.34 

1.11 

0.32 

0.16 

gain 

1.64x 

4.51x 

4.73x 

25.45x 

2.14x 

2.36x 


particularly for the uplink case, which further underscores the 
poor performance of the statically-configured LTE relays for 
uplink traffic. Edge rates, which we define as the lowest 5% 
of rates, also improve notably in all cases. Under proportional 
fairness, the much lower rates of median and edge users in 
the static relay case result in large increases in system utility 
when DTDD allocation is used. 

We also find that, despite adding 4.76 dB of SINR on 
average, the gains of the 4-relay over the 2-relay network are 
only slight for the static TDD case, whereas they are much 
more pronounced (linear for the UL and nearly linear for 
the DL mean rates) under DTDD. This result highlights how 
adding relays provides little benefit when the eNB-RN links 
are bottlenecked, which is evidently the norm for LTE-style 
relays in dense mmW networks like we have simulated. 

It is important to note that the results shown for the static 
LTE-relay case represent the average performance over all 
TDD schemes with 8 usable subframes. This does not illustrate 
the point that the network must initially select a single config- 
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Figure 2: CDF of DL and UL SINR for non-outage UEs after 
beamforming with omnidirectional interference 




(b) Uplink 

Figure 3: CDF of flow rates 


uration, which may result in the ratio of uplink to downlink 
subframes and eNB-to-RN and RN-to-UE subframes being 
very poorly matched to the given the network loading. User 
data traffic in cellular networks is highly variable, so by this 
fact alone, the gains of DTDD should be even more significant 
than shown by these simulations. Also since different relays 
may have traffic loads that vary significantly, a single global 
frame configuration will likely not be appropriate for each 
























































































relay and its users. In future work, we intend to demonstrate 
the full benefits of our scheme by developing a network model 
with time-varying channels and traffic. 

Performance and Optimality: From comparing the per¬ 
formance of our suboptimal DTDD algorithm to the optimal 
allocation, which is computed via brute-force search, we find 
that for the 2-RN/10-UE case, DTDD achieves 13% and 7% 
off the optimal DL and UL mean flow rate and is within 
5% of the optimal utility. Our heuristic took 12 iterations 
(of calling the inner optimizer), on-average, and a maximum 
of 19 iterations to complete, as opposed to 329 for brute- 
force (found by numerically generating and counting unique 
permutations). We separately perform these simulations over 
only 4 subframes since the set of unique configurations we 
must search over increases exponentially with the number of 
subframes and would be intractable for 8 subframes. Clearly, 
the number of permutations of N s f subframes for a single 
BS and Nrn relays is upper-bounded by 2 N s f ( N R N + 1) 5 an d 
while the number of unique permutations is much fewer, it still 
demonstrates that a brute-force approach is computationally 
unmanageable. 

VI. Conclusions 

Multi-hop relaying provides an attractive technology for 
extending coverage in mmW systems which are inherently 
limited by range and blocking. In conventional cellular sys¬ 
tems, the full benefits of relays are difficult to realize due to 
the need to have global, common UL-DL time frame assign¬ 
ments with reserved subframes for relays. However, due to 
directional isolation with high-gain antennas, we have argued 
that duplexing schedules in mmW systems can be selected 
individually on each link. To exploit this flexibility, we have 
proposed a dynamic duplexing TDD system, where duplex 
schedules can adapt to local channel and traffic conditions. 
We have formulated the selection of the duplex schedules as a 
joint optimization with the bandwidth allocation to maximize 
a global utility. The resulting optimization is non-convex, 
but we have developed a computationally efficient subopti¬ 
mal algorithm based on recursive searching within subtrees. 
Simulations of realistic deployments demonstrate considerable 
gains over static allocations, even under fairly uniform traffic. 
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